Anomalies in Universal Intensity Scaling in 
Ultra-relativistic Laser-Plasma Interactions 
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Laser light incident on targets at intensities such that the electron dynamics is ultra- 
relativistic, gives rise to a harmonic power spectrum extending to high orders and char- 
acterized by a relatively slow decay with the harmonic number m that follows a power law 
dependence, m~ p . Relativistic similarity theory predicts a universal value for p = 8/3 up 
to some cut-off m — m* . Results presented in this Letter suggest that under conditions in 
which plasma effects contribute to the emission spectrum, the extent of this contribution 
may invalidate the concept of universal decay. We report a decay with harmonic number 
in the ultra-relativistic range characterised by an index 5/3 < p < 7/3, significantly weaker 
than that predicted by the similarity model. 
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Since the earliest observations of harmonics of laser 
light reflected from target plasmas by Burnett et al. [1] 
and Carman, Rhodes and Benjamin [2], the topic of har- 
monic generation has been one of enduring interest. The 
range of harmonics excited is governed in large part by 
the level of I)? L where / is the intensity of the incident 
light and its wavelength. Levels of I)? L attained in 
recent experiments are such that the laser-plasma in- 
teraction physics extends far into the ultra-relativistic 
(UR) regime. Harmonic orders ranging from ~ 850 for 
peak intensities a little over J 20 = 10 20 Wcm~ 2 to ~ 3200 
when intensities are an order of magnitude greater have 
been reported by Dromey et al. Ref. [3]. In this work 
the decay of intensity with harmonic order is found to 
behave as I m ~ m~ v where 2.2 < p < 2.7 at intensities 
I ~ I20 and 2.5 < p < 3.0 at intensities in the higher 
range. Harmonic decay with p in the range observed is 
consistent with an intriguing prediction by Baeva, Gor- 
dienko and Pukhov (BGP) [4] based on a similarity anal- 

l 

ysis which led them to conclude that the harmonic decay 
from ultra-relativistic laser-plasma interactions is inde- 
pendent of the details of the interaction physics and that 
the value of the decay index p = 8/3 is in fact universal. 

In this Letter we report results from a set of particle- 
in-cell (PIC) simulations in the UR range which suggest 
that the decay of harmonic intensities may be governed 
by plasma effects. In these circumstances we present 
evidence of spectral decay more complex than the uni- 
versal decay proposed by BGP, no longer characterised 
by a unique index p = 8/3 but by indices spanning a 
range 5/3 <p < 7/3. 

The effects that appear in the spectrum reflect de- 



tails of the interaction physics, in particular the plasma 
oscillations generated by bursts of highly relativistic elec- 
trons. We have shown elsewhere that with p— polarised 
light, collective effects not only give rise to emission at 
the plasma frequency and its harmonics [5] but can serve 
as a source of a modulation at the plasma frequency that 
appears in the spectrum for suitably chosen parameters 
[6]. Other work dealing with the effects of plasma oscil- 
lations driven by energetic electrons on emission spectra 
has been reported by Teubner et al. [7] and by Eidmann 
et al. [8]. 

In the simulations carried out to generate the har- 
monic spectra we treat the interaction of a laser pulse of 
length t p with a plasma slab of initially uniform den- 
sity, apart from an interface at the front end. This 
vacuum-plasma interface is characterised by a density 
profile with scale length A that is some prescribed frac- 
tion of a wavelength across the ramp. We have used 
a 1-1 /2D fully relativistic and explicit PIC code, that 
embeds the Bourdier technique [9] to allow for oblique 
incidence. Two vacuum gaps extend from both the front 
of the ramp and the planar rear surface of the plasma 
to the walls of the simulation box to allow for particle 
and wave propagation. When external fields reach the 
boundary and propagate outside the system they are not 
reflected and thus no longer considered in the interaction 
physics. Plasma particles are reflected at the boundaries 
by reversing their velocities. The plasma filled a simula- 
tion box extending over 4 laser wavelengths with 2 x 10 6 
particles distributed across 10 4 grids, giving a ratio of 
grid size to wavelength of ~ 10 -4 . The number of par- 
ticles employed and the choice of the simulation param- 



eters is sufficient to resolve both the Debye length and 
the highest frequency mode with acceptable resolution. 

The initial electron temperature was chosen to be 
100 eV and we used gaussian pulses of variable length 
and profile. Ions were treated as a neutralizing immobile 
background, an acceptable approximation given the fem- 
tosecond time scales governing these simulations where a 
17 fs gaussian pulse is used throughout. The normalized 
quiver momentum ao ~ 8.5 (ho^t) 1 ^ 2 lay in the range 
5-20 with slab densities between 20 n c and 200 n c , where 
n c is the critical density. The reflected emission spectra 
showing the relative strength of the different oscillation 
modes are determined from a time-resolved Fourier de- 
composition of the harmonic contents of the reflected 
electric field at the vacuum gap, normalized to the fun- 
damental. 
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Figure 1: Harmonic spectrum generated from a laser pulse 
incident normally (ab) and obliquely (cd) on a plasma tar- 
get: ao = 10, n e /n c — 40, t p — 17 fs and \l = 1054 nm. 
In (a) the line represents spectral decay P m ~ m~ p with 
a decay index p = 8/3, (b) relative intensity of harmon- 
ics in (a) scaled to m = 100, i.e. P m /Pwo = (m/100) _p . 
The lines correspond to p — 2.1 (dotted), p = 8/3 (solid) 
and p = 3.0 (dashed); (c) harmonic spectrum generated by 
obliquely incident laser pulse showing the effects of spec- 
tral modulation; the line denotes a spectral decay p = 2; 
(d) relative intensity of harmonics in (c) scaled to m = 100; 
the lines correspond to p = 5/3 (dotted), p = 8/3 (solid) 
and p = 10/3 (dashed). 

In Fig. 1 we contrast PIC spectra from runs at UR 
intensities with light incident both normally (Fig. lab) 
and obliquely (Fig. led) on the surface of the plasma 



slab. The plasma density is chosen to be initially uni- 
form with n e = 40 n c and the peak intensity of the 
incident pulse such that ao = 10. Fig. la shows har- 
monic power levels up to m = 1000; harmonics up to 
m ~ 200 are characterised by conversion efficiencies 
> 10~ 6 . The decay in peak harmonic power levels with 
harmonic number m is well represented across the range 
by the index p = 8/3, in agreement with both the simi- 
larity analysis and the PIC spectrum due to Baeva, Gor- 
dienko and Pukhov [4]. At the same time it appears that 
across the harmonic range in Fig. la there is no marked 
break in the spectrum predicted by their model at a crit- 
ical harmonic m* = V8 a jg where y s is the maximum 
value of the relativistic factor associated with the mo- 
tion of the plasma surface and a, the second derivative 
of the surface velocity, is of order 1. For the parameters 
in Fig. 1, m* ~ 40 — 80, adopting the BGP estimate for 

Is- 

Fig. lb presents data from the spectrum in Fig. la 
in which a set of data points is represented over the 
restricted harmonic range m = 40 to m = 400 to af- 
ford a close-up of the macroscopic intensity scaling ap- 
plied across the full harmonic range in Fig. la. We ar- 
bitrarily normalize the harmonic power levels to Pioo 
and find from Fig. lb that across the harmonic range 
(40 - 400) the scaling with harmonic number falls in a 
range 2.1 < p < 3.0. In passing we note that this result 
is not greatly at variance with the scaling reported in 
[3] where 2.2 < p < 2.7 and 2.5 < p < 3.0 (the latter 
applied across a much higher harmonic range), depend- 
ing on intensity. We draw no comparison between our 
simulated scaling and those from Ref. [3] in view of dif- 
ferences in parameters such as the pulse length (17 fs in 
the simulations as against 500 fs in the experiments) and 
without knowing the experimental plasma electron den- 
sity. That notwithstanding, the results for s-polarised 
light in Figs. la,b are broadly consonant with the BGP 
scaling p = 8/3. 

Figs, lc, Id highlight changes in the pattern of har- 
monic emission when p-polarised light is incident at 23° 
to target normal, these spectra showing distinct differ- 
ences from those at normal incidence or for s-polarised 
light. It is at once apparent that on a macroscopic level, 
the spectrum is no longer well characterized by a uni- 
versal decay index p = 8/3. In contrast to the spectrum 
in Fig. la there is now a more or less well-defined break 
in the spectrum in Fig. lc, centred on m ~ 100. We 
emphasize that this feature is quite distinct from the 
spectral break deriving from a validity condition in the 
BGP self-similar model. The break in Fig. lc occurs 
with the onset of modulation; below this the spectrum 



is enhanced over that for s-polarised light (cf. Fig. la) 
by plasma emission with the consequence that the decay 
coefficient across this region is reduced below the BGP 
"universal index" . We show in Fig. Id an intensity scal- 
ing up to m = 400 corresponding to that in Fig. lb, 
with harmonic power levels normalised to -Piooj charac- 
teristic of the distinctive modulation in the spectrum in 
Fig. lc. Below m = 100 decay is less rapid than the cor- 
responding universal decay, with p ~ 5/3, while above 
this the best fit for this choice of parameters is given by 
p ~ 10/3. 
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Figure 2: Modulation-induced anomalies in the univer- 
sal scaling of harmonic spectra from ultra-relativistic laser- 
plasma interactions. The laser pulse with pulse-length 
t p = 17 fs, is incident obliquely on a plane target at an 
angle of incidence 8 = 23° and the wavelength Xl = 1054 
nm. (a) ciq = 5,n e /n c = 20, (b) ao = 20, n e /n c = 40, (c) 
ao = 10, n e /n c — 50, (d) ao = 20, n e /n c — 100. The lines 
indicate a scaling P rn = m~ p , where p = 7/3 in (a), p = 5/3 
in (b),(c) and p = 2 in (d). 

In Fig. 2 we present a range of spectra for additional 
combinations of ao and n e /n c , namely (ao,n e /n c ) = 
(5,20), (20,40), (10,50) and (20,100). All four com- 
binations show spectra with distinctive structure across 
the range and serve to confirm the deviation from BGP 
universal scaling already apparent from Fig. lc. We find 
some variation in the onset of modulation depending on 
particular combinations of ao and n e /n c , with thresholds 
ranging from m ~ 20 — 30 (Fig. 2a, d) to m ~ 50 — 70 
(Fig. 2b, c). Not only is there variation in the onset of 
modulation, there are also notable differences in modu- 
lation frequency across the range of parameters consid- 



ered. This contrasts with our findings in Ref. [6] where, 
for a moderately relativistic range, spectra showed fairly 
regular modulation at the plasma frequency u p , with on- 
set of modulation consistently at ~ 1.5 u p . Despite the 
variation in detail in the spectra in Fig. 2, the decay in- 
dex in each case is tolerably well fitted by p in the range 
5/3 < V < 7/3, markedly less than the spectral decay 
observed in simulations at lower ao [6]. 

The physics differentiating irradiation by p-polarized 
from s-polarized light derives from bursts of Brunei- 
excited electrons [10] that drive plasma oscillations to 
highly nonthermal levels close to the plasma surface. 
Given the steep gradients in density at the front of the 
simulation box, the plasma waves excited are no longer 
uniquely longitudinal but couple strongly to the radia- 
tion field, giving rise to emission at the plasma frequency 
in both back and forward directions, with weaker con- 
tributions at higher harmonics of the plasma frequency 
[5]. Boyd and Ondarza-Rovira further identified a ro- 
bust feature in the spectrum that appeared between the 
plasma line and its second harmonic, typically at a fre- 
quency 1.5 u p , with width of the order of the plasma fre- 
quency. This so-called 'combination line' proved in the 
event to be the first cycle of a modulation that was evi- 
dent in the high harmonic spectrum [6]. Across a wide 
range of combinations of (ao, n e /n c ) Boyd and Ondarza- 
Rovira found that the frequency at which harmonic in- 
tensities were modulated was governed by the plasma 
frequency. A Fourier decomposition of the plasma elec- 
tron density showed correspondingly strong modulation 
at the plasma frequency so that harmonics of the laser 
frequency ul such that m > 1.5 (u p /ul) reflect this 
modulation in density, cf. [6]. 

In the UR-regime the laser-plasma interaction physics 
is strongly non-linear. Not only is the plasma line it- 
self broadened, but its harmonics may also be strongly 
driven and likewise broadened. One consequence of the 
non-linearity is some variability in the order of plasma 
harmonics excited for different combinations of intensity 
and electron density. Another effect of emission at har- 
monics of the plasma frequency I lo v is a widening of the 
range of m well beyond the value m ~ 1.5 {up/u£) at 
which the first minimum appeared in the harmonic spec- 
trum at moderately relativistic energies, Ref. 6. Thus 
one might now expect the threshold for onset of modu- 
lation to be shifted to higher frequencies, typically m* ~ 
(1 — 2) I (up/uLj where I is the highest plasma harmonic 
excited, which may contribute to the variation found in 
Fig. 2. 

While Fig. 2 shows departures from the universal 
p = 8/3 scaling proposed by BGP, fine detail in the 



structure of these spectra is seen more readily over sec- 
tions of the full harmonic range. In Fig. 3a we show a 
Fourier representation of the electron density at a point 
inside the plasma boundary for the parameters used in 
Fig. lc. We see from this that the structure in the 
Fourier decomposition of electron density is reflected in 
the power spectrum. Above m ~ 4lo p /ujl ~ 25 is a 
region extending to m ~ 100 over which the Fourier 
amplitudes show only modest decay, a detail mirrored 
in the behaviour of P m found in Fig. lc, where there is 
a clear discontinuity in the decay index in this region. 

As an illustration of the variability referred to ear- 
lier, we show in Fig. 3b a Fourier density representa- 
tion for the parameter combination in Fig. 2d, i.e. ao = 
20, n e /n c = 100, using a log- linear plot to highlight the 
periodicity of the density modulations. This shows a rel- 
atively broad range across which there is comparatively 
little variation in Fourier amplitudes; between m ~ 80 
and m ~ 140 we find modulation at the plasma fre- 
quency. Above this the modulation frequency switches 
to approximately 2uj p . The triplet that appears over 
170 < m < 230 is the source of broadened emission over 
this range in Fig. 2d, typically ~ Qlo p wide though this 
is difficult to discern given the poor resolution in Fig. 
2d. 





Figure 3: (a) Fourier decomposition of the electron density 
corresponding to the parameters used to generate the spec- 
trum shown in Fig. lc; (b) a log-linear representation of 
Fourier amplitudes across the range 20 < m < 350 for the 
harmonic spectrum in Fig. 2d, where ao = 20, n e /n c = 100. 

The range of values of the decay index in the moderately 



relativistic range 0.5 < ao < 2.0 considered in Ref. [6] 
varied from p~ 4 to » ~ 8/3, for parameter choices such 
that typically S = n e /(aon c ) > 20. In the UR regime 
S < 10 for the three ao values used, ao = 5, 10, 20. Fig. 
4 shows the decay index from PIC runs spanning this 
range of parameters. For ao = 5 the decay index was 
found to be p = 7/3, irrespective of density. For higher 
ao this reduces to p = 5/3, or in a few instances, p = 2. 
For each combination considered, the index appears to 
be distinctly weaker than the universal value proposed 
by BGP. In contrasting the values of p appearing in Fig. 

4 with the BGP universal value 8/3 we emphasise that 
nowise do we attribute even quasi-universality to a par- 
ticular value of p, be it p = 5/3 or 7/3. 

That aside, results from this set of simulations sup- 
port our contention that the decay of intensity with har- 
monic number in the ultra-relativistic regime is more 
complex than allowed by the "universal" decay index 
p = 8/3. The central tenet of the self-similar approach 
is that for constant values of the similarity parameter 

5 = n e / (ao n c ), the laser-plasma dynamics remains sim- 
ilar, since 5-similarity corresponds to a multiplicative 
transformation group of the Maxwell- Vlasov equations. 
While this is indeed the case for s-polarised or normally 
incident light, S no longer serves to characterise the dy- 
namical system uniquely when plasma effects contribute 
to the interaction physics. 
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Figure 4: Deviations of p from 8/3 universality across a 
range of S — n e /(aon c ) for ao = 5, 10,20 and a range of 
electron densities. The double solid circles indicate two or 
more data points at S= 1.5, 2, 4, and 5. 

Thus while the combinations of intensity and electron 
density in Figs, lc and 2a correspond to S = 4, and 
those in Figs. 2c, d to S = 5, the interaction physics, at 
least as far as its footprints in the emission spectrum 
is concerned, does not appear to be similar in either 
case. Note in passing that by similar we intend formal 
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S— similarity, as against the conventional usage; in that 
sense all the spectra are similar. Whereas for the com- 
bination in Fig. lc, p = 5/3, for Fig. 2a, p = 7/3; simi- 
larly p = 5/3 for Fig. 2c as opposed to p = 2 in Fig. 2d. 
We conclude that plasma-enhanced emission together 
with plasma-induced modulation can result in decay of 
the harmonic spectrum significantly weaker than that 
predicted by similarity theory. At the highest intensity 
modelled in this work (tto = 20) we have not found any 
p < 5/3; that notwithstanding, we do not rule out yet 
lower values at higher intensities. Indeed under more 
extreme UR conditions we have obtained PIC spectra 
for which p ~ 4/3. 
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